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We develop a theoretical framework for the exploration of quantum mechanical coherent popula- 
tion transfer phenomena, with the ultimate goal of constructing faithful models of devices for classical 
and quantum information processing applications. We begin by outlining a general formalism for 
weak-field quantum optics in the Schrodinger picture, and we include a general phenomenological 
representation of Lindblad decoherence mechanisms. We use this formalism to describe the inter- 
action of a single stationary multilevel atom with one or more propagating classical or quantum 
laser fields, and we describe in detail several manifestations and applications of electromagneti- 
cally induced transparency. In addition to providing a clear description of the nonlinear optical 
characteristics of electromagnetically transparent systems that lead to "ultraslow light," we verify 
that — in principle — a multi-particle atomic or molecular system could be used as either a low 
power optical switch or a quantum phase shifter. However, we demonstrate that the presence of 
significant dephasing effects on the metastable levels destroys the induced transparency. Finally, a 
detailed calculation of the relative quantum phase induced by a system of atoms on a superposition 
of spatially distinct Fock states predicts that a significant quasi-Kerr nonlinearity and a low entropy 
cannot be simultaneously achieved in the presence of arbitrary spontaneous emission rates. Within 
our model, we identify the constraints that need to be met for this system to act as a one-qubit and 
a two-qubit conditional phase gate. 



I. INTRODUCTION 
A. Quantum Information Processing 

With developments in theory and experiment over the 
last few years, there has been a dramatic growth of in- 
terest in the area of quantum information processing 
and communication, flj 13, 01 There are now strong in- 
dicators that this fundamental research field could lead 
to a whole new quantum information technology in the 
future. Working practical quantum cryptosystems exist 
already. It is known that large (many qubit) quan- 
tum computers, if they can be built, would be capable 
of performing certain tasks (such as factoring large com- 
posite integers or searching) much more efficiently than 
conventional classical computers. Rather smaller (tens or 
hundreds of qubits) quantum processors would be able 
to perform quantum simulations unreachable with any 
classical machine and such processors also have the po- 
tential to extend the working distances and applicability 
of quantum communications. 

There are currently numerous possible routes for quan- 
tum computing hardware. JJ] Many of these are based 
on coherent condensed matter systems. While at present 
such systems generally exhibit less coherence than qubits 
based on fundamental entities (such as ions or atoms), 
the goals are to increase this sufficiently for error correc- 
tion techniques to be applicable and to utilize the po- 
tential fabrication advantage to make condensed matter 
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many-qubit processors. Even if this proves to be the way 
forward, it seems certain that there will be a need for 
separate quantum processors to communicate with each 
other in a quantum coherent manner. Coherent electro- 
magnetic fields are likely the best candidates for realizing 
this goal. In addition, some quantum information pro- 
cessing may be performed directly on photon qubits, us- 
ing non- linear || or hnear quantum optical pr^cesses-^ 
So, photon qubits (or other quantum coherent states of 
the electromagnetic field) certainly have a number of im- 
portant uses for quantum computing. However, they play 
center stage when it comes to communication, because 
the best way to send quantum information over large dis- 
tances is certainly using light, either down optical fibers 
or even through free space. 

Given all this, the study of coherent interactions 
between ligh t and matter is an extremely important 
topic. It is hard to see how large scale quantum 

information technology can emerge without the ability 
to easily interconvert traveling photon qubits and sta- 
tionary matter qubits. 1 12',"l3l Such interconversion would 
enable the construction of quantum networks [14( , where 
communication between matter qubit nodes is mediated 
through photons. The ability to perform quantum gates 
(one- and two-qubit) directly on qubits encoded into pho- 
tons is highly desirable for communication and comput- 
ing. In addition, the coherent storage of photon qubits 
would open up new possibilities for quantum information 
processing. 

Another very important application of quantum phe- 
nomena is to achieve tasks in classical (optical) infor- 
mation processing and communication. These applica- 
tions will likely be the first to impact on future informa- 
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tion technology. In this work we study aspects of coher- 
ent population transfer phenomena and discuss how the 
emergent effects may be applicable for producing large 
optical phase shifts and for switching. These effects in 
themselves are useful for conventional optical data pro- 
cessing. However, if in addition they can be demon- 
strated to work on non-classical input states of light, 
then they have the potential to form the basis of one- 
and two-qubit gates. We present a detailed calculation 
of the relative quantum phase induced by a system of 
atoms on a superposition of spatially distinct Fock states 
(a "dual rail" photon qubit) and identify the conditions 
under which this performs as a useful one-qubit phase 
gate. We also consider the qubit limit of one of the con- 
trol fields applied to the system, and thus how a condi- 
tional two-qubit phase gate can be realized. 



B. Electromagnetically Induced Transparency 

The basic interaction between an atom and a photon 
is weak, which is why photons are so good for communi- 
cation purposes. In order to enhance the interaction, one 
can either seek to use many actual atoms — an ensemble — 
or many "images" of a single atom, in effect created by 
the mirrors of a very high-Q cavity. The latter cavity 
QED approach has produced numerous impressive ex- 
perimental and theoretical results over the last few years, 
flEj l but from the perspective of using atoms to manip- 
ulate light for information processing, we focus on the 
former ensemble approach. |l6j 

Photons, or light in general, can interact strongly with 
an ensemble of atoms, leading to a number of very inter- 
esting effects. The basis of many of these is the idea of 
electromagnetically induced transparency (EIT), a spe- 
cific example of coherent population transfer. |l7| Here, 
quantum interference can be used to effectively cancel the 
would-be absorption in an atomic medium, rendering it 
transparent. Ordinarily, a light signal on resonance with 
an atomic transition |1) |2) (|1) being the lower and 
long-lived level) interacts strongly with the atoms. How- 
ever, if a third long-lived level |3) comes into play and a 
control field resonant with |1) |3) is applied, the atoms 
are stimulated into states which cannot absorb, 0, 
and so in effect the atomic medium is transparent to the 
signal. These so-called "dark states" are coherent super- 
positions of |1) and 1 3). We discuss this effect in detail for 
such three-level atoms. As will be seen, the transparency 
occurs only exactly on resonance. A transparency win- 
dow, necessary for practical applications, can be opened 
by use of four-level atoms. Our focus in this work is 
the use of such EIT systems for the optical information 
processing applications. However, other important and 
interesting directions and applications exist, which can 
be followed up in detail elsewhere. 

The intimate link between absorption and dispersion 
also means that EIT media can be used to manipulate 
the group velocity of light pulses. [23| This can be re- 



duced dramatically below the speed of light in vacuum c 
through reducing the power in the control field and in- 
creasing the atomic density of the EIT medium. Various 
experiments demonstrated the potential of this effect, but 
the significant breakthrough came in 1999 when Hau et 
al. |2j] produced a group velocity of 17 m/s for light 
pulses in sodium vapor. Further developments continue 
to emerge. It is something of an oversimplification to 
regard a light pulse as simply slowing down as it prop- 
agates through an EIT medium. The strong interaction 
of the light with the atoms can effectively be diagonal- 
ized through the concept of quasiparticles (well known in 
condensed matter physics). In this case the propagation 
through the EIT system is described by quasiparticles 
called dark-state polaritons "225 , which are a coherent su- 
perposition of photons and spins. These polaritons move 
at a velocity given by the group velocity of their photonic 
component. An EIT system can therefore act as a "delay 
line" for light pulses, an effect which in itself has signifi- 
cant application potential for communication. However, 
the ultimate limit of this effect is a complete slowing, 
leading to quantum memory for photons. l2^ l25l l2q| 
As the (externally controllable) polariton group velocity 
is reduced towards zero, the photonic component of the 
polariton also reduces. In effect, the quantum state of 
the light field is stored in long-lived spin states of atoms 
within the EIT medium. This clearly has significant po- 
tential for quantum communication and information pro- 
cessing applications. 

In a sense the whole scenario can also be reversed, 
with atomic ensembles (as opposed to the photons) be- 
ing the principal information carriers, with interactions 
being controlled through electromagnetic fields. A num- 
ber of novel quantum phenomena, also with potential for 
communication and processing applications, can arise in 
this case. It is possible to make quantum nondemolition 
measurements on the collective spin degree of freedom of 
atomic ensembles |23| using light. The interaction be- 
tween such ensembles and light has also been employed 
to create a level of entanglement between the collective 
spin degrees of freedom of two atomic ensembles , a 
first step towards QIP with collective spins. Small en- 
sembles of atoms can also exhibit analogous effects to 
the well known solid state mesoscopic phenomena such 
as Coulomb blockade. In very small capacitance sub- 
micron devices the strong Coulomb energy makes just a 
few energy levels relevant, and enables the manipulation 
of individual electronic charges. Likewise in small atomic 
ensembles the strong dipole-dipole interactions make just 
a few energy levels relevant, and it is possible to ma- 
nipulate individual excitations of these spin systems and 
exhibit "dipole blockade" phenomena. |29l ISOl l3ll 13^ 

In this work we develop a theoretical framework for 
the exploration of quantum mechanical coherent pop- 
ulation transfer phenomena, with a particular empha- 
sis on the aspects of nonrelativistic quantum electrody- 
namics that play a strong role in potential classical and 
quantum information processing applications. We be- 
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FIG. 1: Discrete model of a traveling-wave quantum electro- 
magnetic field interacting with an ensemble of atoms within 
a lossless resonator. 

gin by outlining a general formalism for quantum op- 
tics in the Schrodinger picture, and we include a gen- 
eral phenomenological representation of Lindblad deco- 
herence mechanisms. We use this formalism to describe 
the interaction of a single stationary multilevel atom with 
one or more propagating classical or quantum laser fields, 
assuming that this interaction is sufficiently weak that 
the population of the ground state is never significantly 
depleted. In addition to providing a clear description 
of the nonlinear optical characteristics of elect romagnet- 
ically transparent systems that lead to "ultraslow light," 
we verify that — in principle — a multi-particle atomic 
or molecular system could be used as either a low power 
optical switch or a quantum phase shifter. However, we 
demonstrate that the presence of significant dephasing 
effects destroys the induced transparency, and that in- 
creasing the number of particles weakly interacting with 
the probe field only reduces the nonlinearity further. 

II. REPRESENTATION OF THE 
ELECTROMAGNETIC FIELD 

We wish to simplify our treatment of the quantum elec- 
tromagnetic field as much as possible without introducing 
approximations which will result in significant disagree- 
ments with experimental results. Therefore, we begin 
our discussion of field quantization with the model of 
a unidirectionally propagating traveling wave shown in 
Fig. nj As a result of the finite round-trip length L of 
the cavity, the electromagnetic field is comprised of a su- 
perposition of discrete longitudinal modes with different 
wavevectors k, angular frequencies = |k|c, and polar- 
ization unit vectors CkA, where A G {1,2}, that satisfy 
the orthogonality condition e^^iye^^x — (5k',k <5a',a- The 
field interacts with an ensemble of atoms located within 
the cavity; since the resonator is lossless, the only mecha- 
nism able to alter the number of photons in a given mode 
is the stimulated emission and absorption of photons in 
that mode. 

In the laboratory, a typical experiment will most likely 
employ a continuously tunable multi-frequency traveling- 
wave field interacting with an atom in free space. Our 
model can represent such an experiment if two conditions 



are satisfied: 

1. The resonator enforces an amplitude and local 
transverse spatial distribution on the model field 
that is identical to that of the experiment. There- 
fore, the exact number of photons in a given mode 
k is not important when modeling a classical field, 
provided that the mean field intensity (approxi- 
mately given by {ny) Hlui^c/V, where {ny) is the 
mean number of modes occupying the mode volume 
V) and mode shape are preserved. 

2. The cavity length L is sufficiently long that the 
spontaneous emission spectrum of the atom is not 
altered. 

In principle, we can model the tunability of the laser field 
by assuming that the positions of the resonator mirrors 
can be microscopically adjusted. 

In the Schrodinger picture, [t^I we can represent the 
quantized transverse electric field operator as a discrete 
sum over the allowed wavevectors k and possible polar- 
izations CkA^ 

E(r) = ^ekAfkC/k(r) e'""- akA 

kA 
kA 

where the dimensionless resonator eigenfunction Uk{r) 
describes the transverse spatial dependence of a field with 
a characteristic spot size Wk, and satisfies both ?7k(0) = 1 
and the volume normalization condition 

j^d\\U^{v)\' = ^wlL. (2) 

Therefore, the coefficient fk = ^/hiJ^Jeomi]^ has the 
dimensions of an electric field. 

We will apply the creation and annihilation operators 
ajj._^ and OkA to the corresponding Fock (number) states 
through the ladder operator equations 

«kAl"kA> = %/nkA + 1 l^kA + 1) (3a) 
flkA l^kA) = \/^lkA |"kA - 1) (3b) 

The Fock state |nkA) is an eigenstate of the unperturbed 
Hamiltonian[73l| 

H = '^huji^ flk^akA, (4) 

kA 

with the eigenvalue rikA^^k, giving the total energy 
stored in photons of mode kA. In fact, we can ensure that 
we have chosen the correct value for the normalization 
constant 5k by integrating the corresponding expectation 
value of the electromagnetic energy density eoE(r) • E(r) 
over the mode volume. Applying Eq. |(2J), we obtain 

eo / d^r {rikxl E(r) • E(r) |?ikA) = ?^kA hujk, (5) 
Jv 
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as expected. 

In this work, we will pass from the quantum regime to 
the classical limit through the coherent state 
defined in terms of the single- mode (i.e., single frequency 
ui and single polarization A) Fock states \n) as 



<t)) = 



n=0 



(6) 



where a{t) is defined in terms of the mean occupation 
number {ny) as 



a{t) = yflnv) > 



(7) 



Since \a{t)) is an eigenstate of the single-mode destruc- 
tion operator a with eigenvalue a(t), we have 



a\a{t)) = a{t)\a{t)) , and 



{a{t)\a 



t _ 



v{t)\a*{t), 



givmg 



{a{t)\ E(r) \a{t)) = £ U{v) e''^-''a{t) 



C.C. 



C.C., 



where the associated classical field amplitude is 



E{r) - 2 



U{v). 



(8a) 
(8b) 



(9) 



(10) 



Therefore, for a given resonator mode containing exactly 
n photons, we can associate the quantity 2£^/n with a 
classical field amplitude at r = 0. In practice, we 
will allow this classical amplitude — and therefore the 
mean photon number (ny) — to vary slowly in time when 
we study the adiabatic interaction of pump-probe pulses 
in our discussions of semiclassical coherent population 
transfer. However, we must be careful when passing from 
the quantum to the classical regime when the number of 
photons in a pulse is relatively small. In Section llll CI we 
show that weak coherent fields that are used to manip- 
ulate a quantum optical nonlinearity can be disrupted, 
in the sense that the field can no longer be represented 
as a coherent state after the interaction has ceased. In 
all cases, the atom -I- field system dynamics can be cor- 
rectly described by tracking the evolution of individual 
energy (i.e., Fock) eigenstates, and then evaluating the 
sum given by Eq. ^ either analytically or numerically. 



III. ELECTROMAGNETICALLY INDUCED 
TRANSPARENCY 

One of the most striking examples of coherent popula- 
tion transfer is electromagnetically induced transparency, 
where the absorption of a weak probe field by an atomic 
or molecular medium is mitigated by a second control 
field, js^ [s^ In general, we will be concerned with sys- 
tems that are weakly pumped, in the sense that we will 
assume that the perturbed atom(s) will remain primarily 
in the ground state. 




|2> 



(a) Semiclassical energy levels 



CO, 



|2, %> 



CO, 



|2, na - 1) 



\hna-l) 
|2, na - 2) 



(b) Quantum energy manifolds 



FIG. 2: Electric dipole interaction between a two-level atom 
and a nearly resonant electromagnetic field. In the semiclas- 
sical view, the atomic energy levels are separated by the en- 
ergy h u>2i , and coupled by a field oscillating at the frequency 
i^a = (^21 + i^a- In the quantum view, the states of the atom 
-I- photons system separate into manifolds coupled internally 
by resonant transitions; the energy difference hfa between the 
two levels of a given resonant manifold is far smaller than the 
difference huja between two adjacent manifolds. 
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A. Quantum Optics of the Two-Level Atom 

The electric dipole interaction between a two-level 
atom and a nearly resonant electromagnetic field is shown 
sche matic ally in Fig. |21 In the semiclassical view of 
Fig. 2(a) the two atomic energy levels are separated by 
the energy huj2i = huj2 — hwi, and coupled by a field os- 
cillating at the angular fr equen cy uja ~ 1^21- However, in 
the quantum view of Fig. 2(b) the outer-product states 
of the atom + photons system separate into manifolds 
coupled internally by resonant transitions; if we de- 
fine the detuning parameter 



(W2 - t^l) : 



then the energy difference fwa between the two levels of a 
given resonant manifold is far smaller than the difference 
huja between two adjacent manifolds. 

Based on Fig. |2(b)| we write the unperturbed Hamil- 
tonian of the system in the Schrodinger picture 37] as 



Ha ~ huji ail + fiLU2 (J22 + S-'- 



(12) 



where and a are respectively the creation and annihila- 
tion operators for photons with energy htOa and polariza- 
tion vector Ea, and try = \i) {j\ = aj^ are the atomic rais- 
ing and lowering operators. Now, given an atom initially 
in the ground state and photons stored in the (lossless) 
resonator, the two eigenstates of the unperturbed Hamil- 
tonian Hq are |1, Ha) and |2, — 1), with the eigenvalues 



{l,na\Ho\l,na) = huji+nahuja, (13) 



and 



(2, Ua - 1| Ho |2, Ua-l)^ hLU2 + («„ - 1) hUa- (14) 

In the long-wavelength approximation, the electric 
dipole/field coupling interaction Hamiltonian for a sin- 
gle atom at r = is given explicitly bv [33Ll33 .l38j 



V = d- E(0) 

= — (dl2 (712 + d21 (721) 



(15) 



In the basis of the unperturbed Schrodinger eigenstates, 
the matrix elements of the interaction Hamiltonian are 

{2,na -1\V\1, Ha) ^ ea -(121 SaVn^, (16) 
Therefore, if we define the effective coupling constant 



da = ■ d21 £a 

and the effective Rabi frequency 



(^21) then in the unperturbed Schrodinger basis 



{ |l,"a) , 

\2,na-l) } 



we obtain the total Hamiltonian 

H = -h 



ni 



(17) 



(18) 



(19) 



(20) 



where we have subtracted Eq = {l,na\HQ\l,na) from 
both diagonal terms. Note that we are using the bare 
eigenstates (rather than dressed states [s^]) of the un- 
perturbed Hamiltonian in the Schrodinger picture; the 
resulting perturbed Hamiltonian agrees with that of the 
corresponding semiclassical system in the interaction pic- 
ture. 

The evolution of the wavcfunction 

|7/.(i)) = ci{t) |l,n,) +C2(t) \2,na - 1) (21) 
is governed by the Schrodinger eauation|37| 

zhj^m))^Hm)), (22) 

which has the formal solution 

\m) = U{t) 1^(0)) , (23) 
where the evolution operator U (t) is given by 



-,-iHt 



cos {^R t) 



'20j 



cos {flu t) + j jff" ^^'^ (^fi 



(24) 



and 



(25) 



Decoherence phenomena significantly complicate the 



evolution of the two-level manifold of Fig. 2(b) For ex- 
ample, spontaneous emission by the excited atom can 
scatter a photon into the free-space boundary volume 
(the "environment" ) enclosing the idealized resonator of 
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Fig. Hat a rate|3E 



A 



21 



3 TT eo fi. c'^ 



(26) 



causing a transition from the the central manifold of 



Fig. 2(b) to the lower manifold. Strictly speaking, then, a 
completely general model describing incoherent popula- 
tion transfer phenomena must incorporate multiple mani- 
folds. However, in this work, we are primarily concerned 
with applications of nonrelativistic quantum electrody- 
namics to quantum information processing, particularly 
high-fidelity quantum gate operations. In this case, the 
resonantly-coupled manifold of interest is {|1, 1) , |2, 0)}; 
therefore, spontaneous emission by the excited atom 
causes a transition from the product state |2,0) to |1,0), 
a state lying in the lowest manifold. (In Section^ we 
explicitly assumed that the resonator of Fig. ^ does not 
modify the spontaneous emission spectrum, so we do 
not consider spontaneous emission into the cavity mode.) 
This transition effectively destroys the information car- 
ried by the quantum state of the electromagnetic field, 
and therefore a detailed analytical representation of the 
subsequent evolution of the state of the system is un- 
interesting. In practice, we can recover our resonant- 
coupling approximation by extending our product states 
to append an entry indicating whether a photon with fre- 
quency uja has been emitted into the environment, giving 
us a new single-photon basis, extended from Eq. (fl^ to 



{ |l,na,0), 
|2,na-l,0), 

|l,«a-l,l) }. 



(27) 



The introduction of decoherence into our model will 
prevent us from describing the two-level system state us- 
ing the pure vector given by Eq. H21(l . Therefore, in the 
extended basis of Eq. (|27|) . we introduce the correspond- 
ing density matrix 37J of the atom-photon system as 



Pii P\i Pie 

P21 P22 P2e 
Pel Pe2 Pee 



(28) 



We can then use the corresponding total Hamiltonian 

(29) 



H 



ill 

Va 





and an appropriate set of initial conditions to solve the 
density matrix equations of motion js^ 



^- [H,p{t)]-r[pit)], 



(30) 



where we have incorporated damping through the de- 
coherence operator r[/9(f)]. We adopt the Lindblad 
formHIH of r[p(t)], given by 



r(p) 



2 ^ 



7™ {[pLl,,L,-, 



[Ll,L^p\), (31) 



to preserve both positive probabilities and a positive 
semidefinite density operator. The Lindblad operator 
represents a general dissipative process occurring at the 
rate 7™. For example, we can describe the depopulation 
of the product state |j, . . .) (due either to spontaneous 
emission from the atomic state \j) to the environment 
|e), or to photon transmission or scattering losses in the 
cavity) at the rate 7^ using the lowering operator 



L'j = <Te 



and pure dephasing of the state \j, 
using the operator 



L'' 



1 



(32) 

at the rate 7" 
(33) 



where / is the identity matrix. For example, in the case 
j = 2, we have 



000 
000 
1 



(34) 



and 



T" 
^2 



V2 



(ctII — (722 + fas) 



1 

V2 



1 
0-10 
1 



(35) 



Since we have assumed in Section ^ that the resonator 
of Fig. n does not modify the spontaneous emission spec- 
trum of the atom, we have 72 — A21. Therefore, if we 
assume that atomic level |1) is metastable and that the 
cavity is lossless, then we obtain the decoherence opera- 
tor 



r(p) 



721 P12 lel Pie 
721 P2I 722 P22 7e2 P2e 
7el Pel 7e2 Pe2 -"722 P22 



where 



721 — 7l2 



\i2 



12, 



lel 
7e2 



722 
^ 7le 

■ 72e 



^ l'2, 

= 7", and 
1 



272 + 72 • 



(36) 



(37a) 

(37b) 
(37c) 

(37d) 



We substitute Eq. igHl), Eq. and Eq. ^ into 

Eq. ijSni) to obtain 

pn(t) = -2Im[p2i(t)f^:], (38a) 
P22{t) = -l22P22{t) + 2liJi[p2i{t)ni], and (38b) 

P2l{t) = i{Va+ «72l) P2l{t) 

+ina[pil{t) - P22{t)] (38c) 

and then solve for the elements /9ii(i), P22{t), and P2i{t) 
using a "bootstrapping" method for the initial condition 
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Py (0) = Sil Sji. If we assume that the interaction is un- 
saturated (i.e., \fla\ /721 ^ 1) so that pii{t) » P22{t) for 
ah t, and that — to zeroth order in |r2a| — Pii(^) varies 
slowly compared to vJ^o+tIii then P2i(t) will adiabat- 
ically follow pii{t) and assume some quasi-steady-state 
value iD2i- Therefore, if we ignore short-term transients 
we must obtain 



pii(t) ^exp[-2Im(p2if^:)i], 
and subsequently from Eq. H38c|l 

P2i{t)= [l-e'(-"+»^-)*l p2ipn{t), 



where 



P21 



«721 



«721 



-7li 



(39) 
(40) 

(41) 



Collecting results and solving for P22{t), we quickly find 



Pii(i) 

P22{t) 
P2l{t) 



722 T-a 

t^a + 1721 



where 



1 - OT O*^ 2721 

— = 2Im(p21 ^a) 



(42a) 
(42b) 

(42c) 
(43) 



and pee{t) = I- Pii(t)- p22{t). Therefore, Tr {p^) = 1 - 
0(|ria/72ip), and the quasi-steady-state density matrix 
describes a pure state only in the weak- field limit. 

In the limit where the photon number Ua is sufficiently 
small that 1^2^/721 1^ ^ 1, P21 is a valid approximation 
for P2i{t) only when the time t satisfies I/721 ^Ta- 
li is worth noting that this steady state does not result if 
the system parameters are chosen so that the constraint 
72iT'q 1 is not strictly satisfied. In Fig. |31 we have 
plotted the real and imaginary parts of the off-diagonal 
density matrix element P2i{t) as a function of time for 
the case where Va/l2i = 3. Even though 72iTa = 55.6, 
the magnitude of P2i{i) decays significantly on a time 
scale that is only a few times the transient lifetime. The 
characteristic time Tq can be lengthened by increasing 
the detuning Va or by decreasing the Rabi frequency 

The expectation value of the microscopic polarization 
of the atom at r = is simply 



(d(0)) ^ Tr {pd) ^ /52idi2 4- Pi2d2i, 



(44) 



where the first term arises from the annihilation pro- 
cess and the second from the creation process. We 
define the corresponding complex linear susceptibility 
X^^^ {—ujcn'^a) in terms of the macroscopic polarization 
P(0) and the associated classical field amplitude Ea = 
'^^a\/n^ using the expression 



P(0) 



(d(0)) 
V 



C.C.. 



(45) 



0.00 
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(a) Real part 
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^21^ 

(b) Imaginary part 



FIG. 3: Real and imaginary parts of the ofT-diagonal density 
matrix element p2i(t). After transient contributions with fre- 
quency Va and lifetime I/721, both terms adiabatically follow 
pi\{t) and decay with the characteristic time 1/ra. 



where V is the effective mode volume given by Eq. ||2Jl. 
Therefore, we can calculate the complex susceptibility 
using the definition 



> (-Wa, UJa) = ^ , P2I = , (46) 



which, after applying Eq. lO and Eq. to Eq. ||3TJ, 

gives 



Va - il21 2 \ila\ 



Vi + 721 "^oJla 



(47) 



By convention, the real refractive index T]{uja) and the 
linear absorption coefficient K{uJa) are defined in terms of 
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(a) Refractive index: Re x(— (^a, i^a) 



3, 




(b) Absorption coefficient: Im x(— '^'a , i^a) 



FIG. 4: Refractive index and linear absorption coefficient for 
tlie two-level atom shown in Fig. |21 Both the dispersion and 
absorption are largest at = 0. 



the real and imaginary parts of the susceptibihty as 

rf{LOa) = l + Re\x^^'> {-uja,uja)\, and (48a) 



K (uJa) 



V (^a) C 



Im 



x'^'H~^a,u;a)\ ■ (48b) 



We have plotted the real and imaginary parts of Eq. H47|l 
in Fig.^ with units chosen so that Kiiy'a = 0) = 1. Note 
that the frequency with the largest dispersion (at Va = 0) 
corresponds to the frequency with the greatest absorp- 
tion. 

We can increase the susceptibility defined by Eq. H4t)|) 
by increasing the number of atoms placed in the interac- 
tion region shown in Fig.^ We assume that N atoms are 
at rest near r = in a volume that is small compared to 
A zn, where A = ^ttwq is the effective area of the funda- 
mental laser mode at the beam waist, and zr = ttwq/X 
is the Rayleigh length of the mode. Then, in a traveling- 



wave cavity, the Rabi frequency of each atom has a neg- 
ligible spatial dependence, and is given by Eq. lfTH|) . In 
the unsaturated, weak-field cas e, the lower level of the 
quantum manifold shown in Fig. 2(b) has become iV-fold 
degenerate, since any one of the N atoms can be excited 
to the upper atomic level via resonant excitation. We 
therefore define the unperturbed iV-atom basis 



{ |{l},"a,0), 

{2W},n,-l,0 

|{1},71,-1,1) } 



(49) 



where the entry {1} represents the A^-element string 
1, 1, . . . , 1, describing all atoms in the ground state |1), 
and {2^} represents the same string, with the element 
at position k replaced by a '2', indicating that atom k 
has been excited to the upper level |2). In this basis, we 
write the (A^ -I- 2) x (N + 2) density matrix as 



^{1}{1} P{1}{2(1)} P{1}{2(2)} 

P{2(1)}{1} P{2(1)}{2(1)} P{2(1)}{2(2)} 

P{2(2)}{1} P{2(2)}{2(1)} P{2(2)}{2(2)} 

P{e}{l} /5{e}{2(i)} P{e}{2(2)} 



P{l}{e} 
P{2W}{e} 
P{2(2)}{e} 



P{e}{e} 

(50) 

and, if we neglect any interactions between the atoms, 
the Hamiltonian as 



H 



■ ni ni 

fla Va 
Via Va 











(51) 



We can construct the A^-atom Lindblad decoherence 
operator T{p) by allowing each atom to scatter an ab- 
sorbed photon to the environment at the rate 72 (inde- 
pendent of position), and we assume that the pure de- 
phasing of the state | j, . . .) occurs at the rate 7". If we 
repeat the single-atom quasi-steady-state approach that 
led to Eqs. then in the limit IOal/721 <C 1 in the 

noninteracting A^-atom case, we obtain 

Pii(i) — > P{i}{i}(i), (52a) 
P22{t) — > P{2W}{2W}{t)^ and (52b) 
P2i{t) — > P{2W}{i]{t). (52c) 

where Eqs. (|52|) have the same form as Eqs. (|42|) — in- 
dependent of fc — but with the new time constant 

1 _2N^2l\^a? 



ra{N) 



■7I1 



(53) 



Note that we have implicitly assumed that the popula- 
tion of the system ground state has not been significantly 
depleted (i.e., N p22{t) <C pii{t)). 
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Extending Eq. 1)44(1 to the iV-atom ensemble near 
r = 0, we obtain an expectation value of the microscopic 
polarization given by 



N 



(dAr(0)> = Tr (pd) ^ (P{2(*^)}{i}di2 + P{i}{2(*^)}d2i) 
fe=i 

= P21 di2 +pi2 d21, 

(54) 

where, in the quasi-steady-state regime where the time t 
satisfies I/721 < t <C Ta{N), 



~{N} 
P2I 



«721 



«721 



K + 721 



(55) 



Note that ^21^^ scales linearly with N only in the weak- 
field limit where the probability that more than one atom 
has been excited to the upper level is negligible. 

In Section IIV Al we will calculate both the quantum 
phase shift and the scattering rate encountered by a sin- 
gle photon optically coupled with one or more two-level 
atoms in the interaction region of the model resonator 
shown in Fig.^ We can anticipate those results now us- 
ing the definition of the complex susceptibility given by 
Eq. ((46|l in a semiclassical calculation. Suppose that our 
resonator encloses a coherent state with a{t) — y^n^e"*'^* 
and contains atoms in an interaction region of length 
I having the efi^ective volume ^ttWqI. Then we can ap- 
proximate the total phase shift per round trip of the cor- 
responding electromagnetic field as 



A(f : 



Ar] (wa) ujj 



(56) 



where Arj (uja) = rj (uja) — 1 <C 1, and the fraction of the 
circulating electromagnetic power absorbed per round 
trip as 



Pa 



K (uJa) I- 



(57) 



If we substitute Eq. ((55|l into Eq. I|46(l . then we ob- 
tain the TV-atom susceptibility x^^^'f^^ (— ^a), and 
the corresponding refractive index change Arj (cua) — 
iRe [x'l'^^^ (-Wa,Wa)]- li I <S: L, then the round-trip 
time is given by At = L/c, and we obtain for the cor- 
responding classical electromagnetic field amplitude at 
r = 

£;(i) = £:(0)e-*('^»-^")* (58) 
where the complex frequency shift Wa is given by 



Wa 



-'^^XWW J, (59) 



At 2PaAt 2 
or, after applying Eq. ll^ . 



Wa = p^f > ^. 



(60) 



Therefore, in the classical linear-optical limit, the mean 
phase shift accumulated by the ria-photon coherent state 
at frequency Wo after an elapsed time t is 



ip{t) = Rc(Wa)t = -- 



N\na\ 



-ill 



(61) 



and the mean rate at which photons are absorbed and 
scattered by the atoms is, as expected. 



2Im(iy,)na = 2Im (p^f^^!: 



1 



ra{Ny 



(62) 



where Ta{N) is given by Eq. We see that — in the 

limit 722 Ta (N) 3> 1 — this scattering rate is equivalent 
to 722 ^22^^ expected from Eq. H38b|) . This semiclassi- 
cal result is entirely consistent with that of the TV-atom 
quantum calculation found in Eq. (|52|) . as predicted by 
extending Eq. I|38a|) to the A^-atom case. 



B. Transparency of the Three-Level A Atom 



In Section lill Al we developed a formalism to describe 
the quantum optics of a two-level atom. In this formal- 
ism, the effective Rabi frequency Vlj represents the anni- 
hilation of a photon in mode j, and Vj is the diagonal el- 
ement of H corresponding to the (positively) detuned in- 
teraction between mode j and the adjacent atomic tran- 
sition in a given manifold of the energy-level diagram 
of the total system. In this section, we will follow the 
same formalism to analyze the quantum optical proper- 
ties of the three-level atom shown in Fig. |S1 Note that 
the upper atomic energy level |2) and the new metastable 
level 1 3) are coupled by a control field with angular fre- 
quency ujb ~ LU2 ~ 1^3- It is the destructive quantum 
interference established by this control field that results 
in transparency (i.e., vanishing absorption) for the probe 
field at LOa- 

As in the case of the two- level atom, we reduce the 
semiclassical atomic system depi cted i n Fig. 5(a) to the 
quantum manifold given by Fig. |5(b)] We \dF work in 



a manifold corresponding to an atom in state with 
Ua photons in mode a and ni, photons in mode b. We 
again extend our basis to include energy dissipation to 
the environment by appending an entry to each prod- 
uct state, indicating the occurrence of scattering of a 
photon of frequency uJa or LOb- Therefore, the environ- 
ment can be represented by the nonresonant submanifold 
{|1, Tia — 1, TT-b, 1) , |1, 71a, 71b — 1, 1)} that capturcs dissi- 
pated energy and preserves the trace of the density ma- 
trix. Referencing Section IlII Al by inspection in the ex- 
tended unperturbed Schrodinger basis 



{ |1, 71a, Tib, 0) , 

|2,71a - l,71b,0) , 

|3, 71a - 1,71b + 1,0) , 

{|l,nQ - l,nb, 1) , |l,7iQ,7ib - 1, 1)} } 



(63) 
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(a) Semiclassical energy levels 



|1, n^, nj,) 



\3,na-l,ni,+ 1> 



|2, 1,«^> 



(b) Quantum energy manifold 



FIG. 5: Interaction between a three-level atom and a nearly 
resonant two-frequency electromagnetic field. Note that the 
annihilation of a photon of frequency u>k is represented by the 
complex number fifc. 



we obtain the total Hamiltonian 



n* 



Q- 

nt 

ni i^a~i^b 





where we have defined the detuning parameter 

the effective couphng constant 

gb = eb- £b, 
and the effective Rabi frequency 

9b 



vWTT, 



(64) 

(65) 
(66) 
(67) 



and we have subtracted the energy 

£^0 = {l,na,nb,0\Ho \l,na,nb,0) 

from all diagonal terms. 

The dynamics of a strongly-coupled system (e.g., a sys- 
tem having intracavity fields that are sufficiently intense 
that dephasin g ca n be ignored) is often described using 
dressed states, 33] where the submatrix of the Hamilto- 
nian corresponding to the vectors \2,na — l,n;,,0) and 
|3, rta ~ 1, + 1, 0) is diagonalized in the case of per- 
fectly resonant tuning. For our purposes, weak fields and 
linearized (Lindblad) decoherence models generally allow 
the unperturbed eigenvectors to serve as a reasonably 
accurate basis set. For example, if we diagonalize the 
Hamiltonian given by Eq. (|64|l with Va — i^b = 0, we find 
the nontrivial eigenvalues 



rtn = 0, 

51_ = — ri_R, and 

n+ = +nR, 



(68a) 
(68b) 
(68c) 



where ftj^ = \rta\^ + \^b\^ , and the nonzero eigenvectors 



|0) = - 



1(^1 ^ol 



1+) 



V2 



and 



(69a) 
(69b) 
(69c) 



in the basis of Eq. (|63|l . Therefore, if we assume that the 
system is entirely in the ground state at t = 0, we find 



o* 



(70) 



Therefore, in the strongly-coupled case, the population of 
the state |2, Ua — 1, rib, 0) is given by |r2a/^^_R|^ sin^ {^Rt). 
We anticipate, then, that in a weakly coupled system 
with appreciable decoherence this population will remain 
small for all t if |57a/rib|^ ^ 1 and/or \il.R\ ^ 721- 

We now generalize the phenomenological discussion of 
decoherence presented in Section llH Al to describe more 
complex atomic energy-level schemes. We define 7^ as the 
total free-space depopulation rate to the environment of 
product state |j, . . .) arising both from spontaneous emis- 
sion from |j) to all lower atomic levels and from trans- 
mission and scattering losses in the cavity, and 7" as the 
pure dephasing rate for the state \j, . . .). In general, the 
decoherence coefficient jij of the term jij pij appearing 
in the Lindblad decoherence operator given by Eq. 1^1^ 
can be quickly written down using a straightforward set 
of rules: 

1. In all cases, 7^ = jji. 
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2. If z = e and j ^ e, then 



1 



3. If i ^ e and i ~ j, then 



7jj — 7j 



4. If z, J 7^ e and i ^ j, then 



1 



(71) 



(72) 



(73) 



5. The value of the term corresponding to z = j = e 
in r(p) is chosen to ensure that Tr[r(p)] — 0. 

In the following analysis, we assume that atomic 
level 1 3) is metastable so that 73 = 0. We substi- 
tute Eq. linU) and Eq. lEH) into Eq. and seek the 
quasi-steady-state solution in the unsaturated weak-field 
limit |fia/72iP ^ 1, assuming that P22{t) ^ Pii{t) and 
^aP32{t) <C rijp3i(t) for aU t > 0. Then we obtain 
Pii = 1, Pk2 = Plk - 0' Pk3 = Plk - (where k G {2, 3}), 
and 



/O31 



{Va -l'b+ 1731) 



(i^a + il2l){Va - I'h + «73l) - \^b\ 
(i^a + il2l){l^a - Vb + «73l) " |f^f)|^ ' 



r, (74a) 



(74b) 



with pi2 = P21 and P13 = pt^i, where 721 is given by 
Eq. l|T7a|) and 



731 = 7i + 73 • 



(75) 



Since Eq. (|38a|l remains valid for the three-level atom -t- 
photons system, we can follow the same bootstrap proce- 
dure to obtain the approximate solution for pii{t) given 
by Eq. and 



P2i{t) = P21 (1 
P3i{t) = P31 (1 



-721*^ 



pii{t), and 
Pii(i)- 



(76a) 
(76b) 



Therefore, the steady-state solutions given by Eqs. l(7^ 
are valid at any time t where the laser parameters have 
been chosen to allow the inequality 



1/721,1/731 « t « Ta = [2Im(p2lf^:)]"' 



(77) 



to be satisfied. 

Substituting Eq. (|74a|l into Eq. we obtain the 

susceptibility (linear in Ea) 



{Va -l^b + i73l) 



(78) 



(Va + ij2l){l^a - Vb + 1731) - \^bY 



When 731 — > 0, the new detuning terms in Eq. H78|l have 
remarkable implications for both the refractive index and 
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(a) Refractive index: Re x(— i^a, '.''a) 
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(b) Absorption coefficient: Imxi—i^a, i^a) 



FIG. 6; Refractive index and linear absorption coefficient 
for the three-level atom shown in Fig. |3 with vt — jai — 0. 
Both the dispersion and transmission window are sharpest 
when |f2ft/72i|^ < 1. 



the absorption coefficient given by Eq. H48|l . In Fig. 
we assume that the control field has been tuned to the 
resonance frequency of the |3) |2) atomic transition so 
that = 0, and we plot the real and imaginary parts of 
the complex susceptibility as a function of the normal- 
ized detuning 1^0/721 for several values of the normalized 
control Rabi frequency |f2f,|/72i. When = 0, the 
absorption and dispersion curves reduce to those of the 
two- level atom shown in Fig.^ However, when \ > 0, 
near Va — the absorption vanishes completely over 
a frequency range with a full width at half-maximum 
(FWHM) of 



-7|i 



721- 



(79) 



Over the same frequency range, the dispersion of the re- 
fractive index drj{uj)/duj is proportional to 72i/|f2{,p, re- 
sulting in a significant increase in the group refractive 
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index 77g(w) = ri{uj) + lu dri{ui) / dkv and a corresponding 
reduction in the group velocity c/r]g{uia) at frequency uja- 
However, the dispersion cannot be made arbitrarily large 
by reducing the amplitude of the control field, because 
in the limit |ri6/72ip ^ 1 the width of the transparency 
window given by Eq. ^ becomes Ai^a = 2\nb\'^ /-f2i- 
Instead, the magnitude of the control field must be cho- 
sen to allow all significant spectral components of the 
probe pulse with carrier frequency LUa to be transmitted 
with the maximum possible dispersion. (In an inhomo- 
geneously broadened medium such as a dilute gas, only 
those atoms with Doppler-shifted resonance frequencies 
that are coincident with uJa influence the group velocity 
of the probe beam. [13, EH) 

Note that the transparency predicted by Eq. H74a|) 
arises whenever z/q = i^b, generating a pathway for a "vir- 
tual transition" between the unperturbed atomic energy 
levels |1) and |3). In the steady-state case, the absorption 
by the atom of a photon of frequency cua never occurs, 
in the sense that a measurement of the state of the sys- 
tem will never find the atom in the level |2). Instead, 
the control field creates a coherent superposition of the 
|1) |2) and |3) |2) paths in Hilbert space such that 
destructive interference effectively reduces the |1) |2) 
transition rate to zero. 

We can estimate the Rabi frequency resulting from a 
particular choice of experimental parameters by noting 
that both the spontaneous emission rate ^21 given by 
Eq. 12()|) and the Rabi frequency given by Eq. H18|) depend 
on the electric dipole matrix element d2i. If we assume 
that ea ■ d2i = |d2i|, then wc find that 



I fin I = 7-^21 AuJrnn, 



(80) 



where cr^ = 3A^/27r is the resonant atomic absorption 
cross section at wavelength Aa = 27rc/a;2i,|33 A = 
^ttwq is the effective laser mode cross-sectional area, 
and Aujr = 2'kc/L is the free spectral range of the ring 
resonator. If we were simulating the adiabatic interac- 
tion of a pulsed laser field with a stationary atom, then 
AuJr would represent the bandwidth of the pulse profile 
function. 34v 42. 43, 44, 45 ] although more complex tran- 
sients can arise when the interaction is non-adiabatic. 
In the weak-control-pulse case where \flb\ /7I1 ^ 1, we 
would require that AcUr < Ai^a — 2 \ flb\ /721 for maxi- 
mum transmission of the signal pulse at frequency coa- If, 
as an example, we choose Acur = \nb\ /721, the Rabi fre- 
quency required to open a sufficiently large transparency 
window is 



A 721 l^al 
0-a ^21 na 



(81) 



Assuming realistic optical focusing parameters, in free 
space the optimum value of the ratio aa/A is about 
20% for Gaussian beams, ^3, ^3 but in waveguides it 
can approach unity. neglect dephasing and set 

li = 72 = 73 = 0, then 721 = ^72 = 5(^21 + A23), 



and 721/^21 « 1 if the branching ratio A23/A21 is ap- 
proximately unity. Therefore, in order to open the trans- 
parency window just wide enough to admit every photon 
in the probe pulse, we must have \flb\ ~ 10|f2a|/Y^n^ 
in free space, and about half that value in a waveguide. 
As an example, for interactions of single probe photons 
having |r^a|/72i ~ 0.03 with rubidium atoms (with a 
lifetime of level |2 = 683/2, F — 3) of 27 ns), the above 
parameters require |f^bP/72i ~ 0.1, giving a FWHM of 
the transparency window of about 1 MHz. Therefore, 
only pulses with durations longer than 300 ns can safely 
propagate through the window. 

The assumption that 731 ^ is generally valid in di- 
lute gases, where decoherence arises primarily from spon- 
taneous emission, although the introduction of Doppler 
broadening alters the frequency dependence of the sus- 
ceptibility given by Eq. (EUl.lialiJ The ultraslow group 
velocity of light propagating through an appropriately 
prepared gas sample ,21] can be reduced dynamically to 
zero, allowing probe pulses at frequency uja to be stored 
in the resulting atomic coherence and subsequently re- 
covered after a controllable time delay. [49] Both trans- 
parency and storage of light has been demonstrated 
in a solid-state system (Pr:YSO) where cooling to hq- 
uid helium temperatures reduces the dephasing interac- 
tions between atomic levels ]1) and ]2).!50] However, in 
semiconductor-based systems, where the possibility that 
731 ^ may be higher given the dependence of the off- 
diagonal decoherence rates on a common set of Lindblad 
parameters, these conclusions become less certain. [5l| 
For example, in units where the absorption of the equiva- 
lent two-level system (corresponding to jJlfc] = 0) is 1, the 
absorption coefficient of the general three-level system at 
i^a = is given by 



K(i^a = 0) cx 



721731 



■ 721731 



(82) 



Hence, if we assume that jrib] /721 — 0.1 and 731/721 = 
0.01, we obtain an absorption coefficient that is 50% of 
the corresponding two-level value, and the system is no 
longer transparent. Similarly, the dispersion of the gen- 
eral three-level system is 



dr]{iy, 



.^ = 0) 721 (|f^b|^-7|i) 



dfa 



(I 



nb 



■ 721731 



(83) 



indicating that we must have \ flb\ > 731 to maintain a 
positive dispersion and a reduced group velocity for the 
probe field. For 731 7^ 0, the optimum value of the control 
Rabi freque ncy yielding t he largest dispersion is given 
by nb = \/72i73i + 27I1 =0.1 when 731/721 = 0.01, 
resulting in a dispersion that is reduced by a factor of 
four relative to the case 731 = 0. 

If we follow the approach of Section IIII Al and extend 
our single-atom three-level analysis to include N atoms, 
we draw essentially the same conclusion as in the two- 
level case: the susceptibility is enhanced by a factor of 
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N in the unsaturated, weak-probe-field limit. Adopting 
the notation of Section llll Al in the unperturbed basis 



{ |{l},7lQ,71b,0) , 

{2(^-)},na-l,nb,0^ 

{3(^-)},7i,-l,nt + l,0), 
{|{l},«a - i,nb,l) , \{l},na,nb - 1,1)} } 
we have the A^-atom Hamiltonian 



H 



(84) 
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(85) 



As in the case of A^ noninteracting two-level atoms, the 
off-diagonal density-matrix elements for each atom are 
given by Eqs. 174(1 . and Eq. ((54(1 holds in the three-level 

case, with p!2^^ = N p2i as before. The aggregate suscep- 
tibility is enhanced by a factor of A^ only if the conditions 

P22^ ^ p\i^ ^^"i 731 ~* satisfied; otherwise, when 
N ^ \ and Va — i^b = 0, the transparency window is 
completely destroyed and the susceptibility reduces to 
that of A^ two-level atoms. However, when these condi- 
tions are met, the absorption coefficient approaches zero 
as Va~Vb 0, and the group refractive index is increased 
by a factor of N . This property allows us to compensate 
for the need to open a sufficiently large window to allow 
a pulse to propagate through the interaction region with- 
out significant loss by increasing the number of atoms in 
that region. 



C. Tunable Transparency of the Four-Level J\f 
Atom 

The capability to transmit a particular probe field with 
very high fidelity and extraordinarily slow group veloc- 
ity is enabled by the inclusion of the second control field 
in Fig. |5| However, since the probe frequency providing 
the greatest transparency results in a differential refrac- 
tive index A?7(a;a) that is precisely zero, this system can- 
not be used to generate a significant relative phase shift. 
More precise deterministic control of the phase and in- 
tensity of the probe pulse can be obtained by adding 
the third field shown in Fig. [71 coupling the metastable 
atomic energy level |3) with the upper energy level |4). 
As in the previous sections, we reduce the semiclassical 
atomic system depicted in Fig. |7(a)| to the quantum man- 
ifold given by Fig. 7(b)| and we will work in a manifold 
corresponding to an atom initially in state 1 1) , with ria 
photons in mode a, rib photons in mode 6, and ric photons 



in mode c. We again extend our basis to include energy 
dissipation to the environment by appending an entry to 
each product state, indicating the occurrence of scatter- 
ing of a photon of frequency Wq, Ub, or ujc- Therefore, 
the environment should be represented by a nonresonant 
sub-manifold that captures dissipated energy and pre- 
serves the trace of the density matrix, resulting in the 
unperturbed Schrodinger basis 



{ |l,na,nh,nc,0) , 
\2,na - l,nh,nc,0) , 
|3, na -l,nb + l,nc,0) , 
|4, ria -l,nb + l,nc - 1,0) , 
{ |l,nQ - l,nb,nc,l) , 
\l,na,nb - 1,71c, 1) , 

\l,na,nb,nc - 1,1) } }. 



(86) 



Referencing Section flll Bl by inspection we then obtain 
the total Hamiltonian 





- 


n* 








0- 








^b 








h 





ni 


Va - Vb 















Vic Va 


-Vb^ 


-Vc 




















H ^ -h ni Va ~ Vb , (87) 

Q,c Va — Vb + Q 

. 0. 

where we have defined the detuning parameter 

Vc = LOc- {i^A - t^a) , (88) 
the effective coupling constant 

gc = ec- di3,£c, (89) 
and the effective Rabi frequency 

= j-^c, (90) 
and subtracted the energy 

= (1,71a, 716, nc,0|i/o |1, 71a, 71b, 71c, 0) 

from all diagonal terms. 

Again, we follow the general decoherence conventions 
established in Section IIII Bl to enumerate contributions 
to the decoherence operator using Eq. (|^ . We de- 
fine 7^ as the total free-space depopulation rate to the 
environment of state |4, ria — 1, 7ib -I- 1, ric ~ 1, 0) and 7" 
as the corresponding pure dephasing rate. We substi- 
tute Eq. (jSTI) and Eq. into Eq. and seek the 
quasi-steady-state solution in the unsaturated weak-field 
limit |^2a/72iP < 1, assuming that pi2{t) <C pii(i) , 
f^a/332(i) < fllpz\{t),^lPAx(t), and flaPA^it) < ^cP-i\{t) 
for all t > 0. We neglect all contributions to density ma- 
trix elements of order \flaf' and higher, and we obtain 

P\\ = 1, Pk2 = P2k - 0> Pk3 = Plk - 0, pki = Plk = 

(where k e {2,3,4}), and 
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(a) Scmiclassical energy levels 



|3,«fl- l,«c> 



|2, l,«^,«c> 



|4, na-l,n^+ l,nc-l> 



(b) Quantum energy manifold 



FIG. 7: Interaction between a four-level M atom and a nearly 
resonant three-frequency electromagnetic field. Note that the 
annihilation of a photon of frequency oj^ is represented by the 
complex number 57^. 



P21 = 



P31 = 





{Va -Vb+ «73l)(j^a -Vb + l^c + »74l) - 






[Va + «72l) 




- [Va 


- Vb + Vc+ «74l) \^b\^ 






{l/a - Vb + Vc+ «74l) 







{Va + 1721) {^a -Vb+ il3l){^a - I^b + t^c + «74l) ^ \^c\^ - {^a - Vb + J^c + ?'74l) \^bf' 

^5 



{Va + «72l) [l^a - l'b + h3l){t^a - Vb + I'c + «74l) - \^c\^ - {^a - Vb + Vc + «74l) l^^bl^ 



(91a) 
and (91b) 
(91c) 



where the remaining ofT-diagonal density matrix elements are pik = p%i [k G {2,3,4}), 721 is given by Eq. (|37a|l . 
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731 is given by Eq. (f75|l . and 741 is defined as 

741 = ^74 + 7^ + 74 (92) 

Since Eq. (|38a|l remains valid for the four-level atom -t- 
photons system, we can follow the same bootstrap proce- 
dure to obtain the approximate solution for pii(t) given 
by Eq. and 



Pkiit)^Pki (l-e-^"*)pii(t), 



(93) 



where fc g {2, 3, 4}. As in the three-level case, the steady- 
state solutions given by Eqs. H91|l are valid at any time 
t where the laser parameters have been chosen to allow 
the inequality 

1/721, 1/731, 1/741 « i « Ta = [2lm{p2ini)]-^ (94) 

to be satisfied. 

We see immediately that Eq. (I91a|l reduces to Eq. H74a|) 
if \/^c + 741 ^ l^cP- Therefore, in this limit the four- 
level system of Fig. [7| remains at least approximately 
transparent. However, in general, the absorption and 
dispersion curves shown in Fig. El can be significantly 
modified by a nonvanishing radiative coupling between 
the atomic states |3) and |4). For example, when all 
coupling fields are resonantly tuned to the corresponding 
transitions (i.e., i^a = = = 0), we have 



P21 



(731741 + |f^c|^) 
721 (^731741 + |f^c|^) + 741 l^^bT 



(95) 



Since the complex part of p2i is related to the absorption 
coefficient by Eq. H48|l. we see that we have lost perfect 
transparency when the frequency coc is resonant, even if 
731 = 0. However, the presence of the fourth atomic 
level and the corresponding nearly resonant electromag- 
netic field with frequency lUc introduces an extremely 
large Kerr-like nonlinearity even at relatively low light 
intensities. JO, 52^ 53, 54, 55, 5(^^^ In particular, we 

can demonstrate that the four-level system of Fig. [7| can 
behave as either a low-energy quantum switch |58i. i59i . 16OI 
Isij or a quantum phase-shifter, [S^ properties which can 
be harnessed for both classical and quantum information 
processing and distribution. [TlL I2II2II29II6I l63l| 

We assume that the magnitude of the Rabi frequency 
fie satisfies the inequahty \flc\^ <C 7I1, so that we can 
define a complex third-order susceptibility for the |1) < — > 
1 2) transition by analogy with Eq. H45|l as 



P(0) = V-i(p2idi2+Pi2d2i) 



(96) 



+ ^a^X*-^^ (-"^a, t^c, -t^c, Wo) \Ec\^ Ea + C.C. 



where x'"'^^ (— Wa,Wa) is the three-level linear susceptibil- 
ity given by Eq. (|78|l . If we expand p2i in a power series 
about 1 57c I ^ — 0, we obtain 



[Va ~ l'b + yc+ «74l) (i^a + i72l)('^a - + 4731) - 



(97) 



As in Section lillBI we set Vb = 731 = 0, we choose the 
values of |ri6|/72i = |f7c|/74i = 0.1 (consistent with weak 
fields), and for simplicity we select 741 — 721. In Fig.|S|we 
plot the real and imaginary parts of the total susceptibil- 
ity as a function of i'all2i for three values of the relative 
detuning Vc/lii- As we discussed above, when — 0, 
the probe field at LOa will be strongly absorbed, and when 
VcIivl ^ 1 transparency is largely restored, but there is 
a nonzero contribution to the refractive index at uja ■ We 
can see these effects more clearly by setting = in 
Eq. H97I) and then varying j/c, as is done in Fig. |21 We 
note that when v,. = 0, the presence of a photon at ojc 
(i.e., Ificl 7^ 0) closes a "quantum switch" that causes 
the absorption of the probe photon, |5^|63,E| and when 
Vc/lAi ~ 30, the atom acts as a "quantum phase shifter" 
that is largely transparent but shifts the relative phase 
of the probe photon, [sll We explore the latter property 



of the four-level system in the next section, where we ex- 
plicitly calculate the corresponding applied phase shift of 
a coherent superposition of single-photon states. 



Following the approach of Section IIIIBI an extension 
of the single-atom four-level model to include N atoms 
in the interaction region reveals that the susceptibility 
is enhanced by a factor of TV only in the unsaturated, 
weak-probe-field limit. Again adopting the notation of 
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Section lill Al in the unperturbed basis 
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(98) 



FIG. 8: Refractive index and linear absorption coefficient 
for the four-level atom shown in Fig. with 731 ^ 0. Both 
the dispersion and transmission window are sharpest when 
l!^b/72i|^ < 1- 



we have the A'^-atom Hamiltonian 
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(99) 



As in the case of A^ noninteracting two- level and three- each atom are given by Eqs. (|91|l . and Eq. (|54|) holds 
level atoms, the off- diagonal density-matrix elements for \^ the four-level case, with p^f ^ = N as before. The 
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FIG. 9: The behavior of the third-order susceptibiUty shown 
in Fig. |H| at = as a function of the normahzed detuning 
Vcl^A\- Note that as z^c ^ 0, a large relative absorption 
arises at the probe frequency LOa (a "quantum switch"), and 
that large values of Vc cause a phase shift that is substantial 
relative to the small absorption (a "quantum phase-shifter"). 



aggregate susceptibility is enhanced by a factor of iV only 

if < Pii a-nd 0. 

Our discussions have emphasized that the resonant 
nonrelativistic quantum electrodynamic interaction of 
the four-level system of Fig. [7| results in the generation 
of a giant third-order optical nonlinearity, commonly de- 
scribed as a Kerr nonlinearity in the literature. [T^ l^l, 
El El El El El m Strictly speaking, an effective 
Kerr Hamiltonian with the form 

Hy^^^^KW a'ac^c (100) 

causes the Fock state |V'(0)) = |"-a,f^c) to evolve as 

= e-^H'at-tct 1^(0)) = e-"»"=^ n,) , (101) 

where (/> = Wt. Therefore, if the evolution of the four- 
level system shown in Fig. [7| exhibits this Kerr behavior, 
then we can claim that the corresponding nonlinearity is 
in fact a Kerr nonlinearity. However, at low light levels, if 
one of the optical transitions is driven by a weak coherent 
state (rather than a Fock state), we can show that the 
structure of this nonlinearity is not strictly of the Kerr 
type unless \^bf 3> |ila|^ , |17c|^- We begin by assum- 
ing that all dephasing rates are zero, and by solving the 
Schrodinger equation for the basis set given by Eq. H8t)|) 
under the influence of spontaneous emission only. Using 
our adiabatic bootstrap approach, we find for the case 
I'a = f fc = that the A''-atom ground state evolves as 



|{1} ,nQ,nb,nc) 



-iWt 



|{1} ,na,n&,nc) , (102) 



where 



W 



+ i ^741 + 721 |l^c|^) 



(103) 



and 7/ji = 7^/2. Note that when the inequality 



721 741 



> 



721 



+ 



I a 



741 



(104) 



is satisfied (equivalent to the assumption i^c/74i ^ 1 in 
the simple case where /721 ~ \^cf /741), the prob- 
ability that a single photon with frequency Wq will be 
scattered by the atom becomes vanishingly small. There- 
fore, the evolution of the atomic ground state during 
a prolonged interaction with the compound Fock state 
jria, nf), Uc) will be governed primarily by the real part of 
Eq. H103|) : since |ria|^ cx Ua and \flc\'^ cx ric, the evolu- 
tion of Eq. (|102|l has the Kerr form of Eq. H101|) with the 
nonlinear coefficient 



W 



A^|»aP|»c|' 



(105) 



where the vacuum Rabi frequencies are given by = 

Let us now replace the nb-photon Fock state with a co- 
herent state parameterized by ah and assess the evolution 
of the corresponding unperturbed ground state 



«r irn \ (106) 

—j= |{l},na,nfc,nc) . 



Since each unperturbed eigenstate evolves according to 
Eq. H102|) . after a time t we find 



\m) 



e 2 



(107) 



where (t>{t) = Wt and 



W 



N\ha\^\h,\ 

j/c|f^fcP|ab| 



(108) 



Note that \'ip{t)) is not a coherent state unless \ab\ ^ 1, 



for which |V'(i)) 



1-0(0)). In Fig.Uni we have 



used Eqs. I|106I) - (|1U7|) to numerically evaluate the inner 
product |(V'(i)|V'(0)) 1^ for several values of the net phase 
shift 0, assuming that ria = 1 and ric — 5. Note that for 
large phase shifts the inner product differs significantly 
from unity even when \ab\^ « 1000; in fact, for = tt, 
|(V'(t)|V'(0))|^ > 0.99 only if > 2.5 x lOl There- 

fore, only when the coupling field driving mode b closely 
approximates a classical field does EIT provide a true 
cross-Kerr nonlinearity. 

Nevertheless, we can appreciate the magnitude of the 
optical Kerr nonlinearity even at low light levels by es- 
timating the parameters included in Eq. I|108l) . Let us 
assume that « |f^c|, and compute the phase shift in- 
duced by a system of 1000 non-interacting atoms in the 
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1000 




FIG. 10: Numerical evaluation of the coherent-state-coupled 
inner product \(ip{t)\ip{0))f for several values of the net phase 
shift assuming that the input Fock state has Ua = 1 and 
He = 5. Note that for large phase shifts the inner product 



differs significantly from unity even when \at\ 



1000. 



FIG. 11: A model Mach-Zehnder interferometer illustrating 
an architecture for a quantum phase-shifter using the four- 
level atom described in Section llVBI 



case where |ab|^ ~ 25. Using Eq. (|5n|l as a guide, we 
estimate aa/A ~ 20%, and we assume a unit branch- 
ing ratio for spontaneous emission from atomic level |2) 
so that 721 ~ ^21- If we let the Fourier-limited pulse 
duration be 2tt/Auj, then after the pulse has interacted 
with the atoms we obtain a phase shift of approximately 
0.1 radians. This shift is about seventeen orders of mag- 
nitude larger than the corresponding value provided by 
a standard Kerr cell. [63, l6a| 



IV. QUANTUM INFORMATION PROCESSING 

We wish to assess the utility of coherent population 
transfer for the creation, transmission, reception, stor- 
age, and processing of quantum information. In partic- 
ular, we must evaluate the time dependence of coherent 
superpositions of discrete states of the atom + photon 
field that are (at least in principle) easily distinguished 
by direct detection of a photon with energy huja- For ex- 
ample, consider a system that is initially in a pure state 
consisting of a superposition of two manifold states, such 
as 



V2 



(|l,0) + |l,na)) 



(109) 



where 



1. 1 1,0) represents the atom in the ground state and 
zero photons in the resonator of Fig. ^ and 

2. |1, Ua) represents the atom in the ground state and 
Ha photons in the resonator. 

If we subsequently apply the unitary phase shift operator 
^(ip) = e*"^"^" to 1-0), then we obtain the result 

1 



V2 



(|l,0) + e»"'^|l,n,)), (110) 



where we note that each Fock photon contributes equally 
to the total accumulated phase. Similarly, if we begin 
with a superposition of an empty resonator and a coher- 
ent state ll.a), we find 



1 



|^')=*(¥')l^) = ^(|l,0) + |l,ae'^)). 



(111) 



in accordance with our intuition for classical fields. A 
simple physical implementation of such a dual-rail co- 
herent superposition could be provided by the Mach- 
Zehnder interferometer shown in Fig. 1111 In one arm of 
the interferometer, the single four-level atom represented 
by Fig. [3 is prepared using \il.c\ > to provide a phase 
shift at the probe frequency uja while remaining largely 
transparent and dispersive. In the second arm, jfid = 0, 
and the system is tuned to match the absorption and 
dispersion provided by the atom in the first arm, allow- 
ing the interferometer to remain time-synchronous. In 
principle, we can simply use the real and complex parts 
of the susceptibilities given by Eq. H78|l and Eq. (|97|l to 
determine the classical absorption and group velocity re- 
duction provided by each system, as is done in standard 
treatments. 0,0,3 In practice, however, we must be care- 
ful to demonstrate that the interaction of either arm with 
a photon at the probe frequency that has entered the in- 
terferometer at the input port will entangle the quantum 
mechanical paths of that photon with each other but not 
with either of the atoms. 

In this section, we solve the density matrix equation 
of motion given by Eq. H30() for product states that in- 
clude the additional kets enumerated above, and we seek 
expressions for density matrix elements that allow us to 
directly read out the phase (p in Eq. (|110|l and Eq. Hlll|) 
in terms of experimentally determined parameters such 
as Rabi frequencies and laser detunings. We will dis- 
cover that constraints must be placed on possible values 
of these parameters in systems suffering from decoher- 
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ence because of the necessity of maintaining either high 
fidehty (or low entropy) in systems without active quan- 
tum error correction, or (equivalently) high data rates in 
corrected systems. 



A. The Two-Level Atom 

Following the example of the dual-rail state introduced 
in Eq. H109() . we wish to further extend the basis of 
Eq. l|77|l to include the possibility that all Ha photons 
have followed another quantum trajectory, and are not 
found in the interaction region containing the two-level 
atom(s). We add one element, and one state vector, 
to include the atomic variables and the environment in 
Eq. fillip : 



{ |l,na,0,0), 
|2,n,- 1,0,0), 
|l,«a-l,l,0), 

11,0,0,71,) } 



(112) 



Now we can rewrite the two-level density matrix given 
by Eq. (^5)) in this basis as 



Pll Pl2 Pie Pio 

P21 P22 P2e P20 

Pel Pe2 Pee PeO 

.POl P02 POe POQ. 



(113) 



We can then use the corresponding total Hamiltonian 



H = -h 



■ 



. 








0" 





0. 



(114) 



to solve for the density matrix element pio{t). We in- 
clude in our model the same phenomenological decoher- 
ence mechanisms introduced in Section IlII Al and we do 
not introduce new dephasing processes between the su- 
perposed loaded and unloaded resonators. 

At t = 0, we assume that the system is in the pure 
state superposition 

|V(0)) = -^(|l,0,0,n„) -Ml,na,0,0)), (115) 

and we wish to identify a later time t = ti (if possible) 
where the system state vector has evolved to the pure 
state 

= ^ 0, 0, na) + e'^(*i) |1, na, 0, 0)) , (116) 

i.e., a state where both the atomic state and the environ- 
ment can be factored out of the Hilbert space, leaving a 
nonzero relative phase difference between the remaining 



basis vectors. At each of these times, the density matrix 
will have the form 



Pit) 



ie*'^^*)' 








-ifW 



(117) 



where Lp{0) — establishes the initial condition. 

In the absence of decoherence, the nonzero density ma- 
trix elements are quickly found to be 



1 



Y- sin^ (nut) 



P22{t) = ^sin2(r!^t), 
P2i{t) = p\2{t) = i-^si-a{VLRt) 

COS (riflt) + i '^ sin (flRt) 



Pw{t) 



P20(t) 



COS {ytf>t) — i ° sin {O.Rt) 

2i 2 r 



-e^''"*sin (nRt), 



(118a) 
(118b) 

(118c) 

(118d) 
(118e) 



where flR is given by Eq. H25|l . It is clear that p{t) has 
the form given by Eq. 1)1171) at the times tq = q-k/VIr, 
where g is a nonnegative integer. At these times, the 
argument of pwit) is given by 



2VLj 



(119) 



where we have chosen the sign of Lpq to be consistent with 
that of the argument of piQ it) when t is small. When the 
system is undamped, a field with the detuning 



2{q-l) 



acquires a relative phase of — tt at the time 



V2^^TT 



(120) 



(121) 



Therefore, for a given value of jJlal applied to an un- 
damped atomic system, a resonantly tuned field with 
q = 1 and Va = obtains a — tt phase shift earliest. 

In the presence of decoherence, we seek the quasi- 
steady-state solution to the density matrix equation of 
motion defined by Eq. H30() . now emphasizing the single 
matrix element pio{t). We build the Lindblad decoher- 
ence matrix operator by applying Eq. H31|l to the new 
density matrix Eq. H113|) . and we extract the system of 
coupled linear differential equations 

Pioit) = -jiQpioit) + inlp2o{t), and (122a) 
P2oit) = iii^a + ij2o)P2oit) + ii^aPw{t), (122b) 
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where the decoherence constants are 

710 = 7": and 
1 , 

720 = -jl2+l2- 



(123a) 
(123b) 



energy levels is small enough that 72 ^72- In this limit, 
dephasing can be neglected, and the argument of the den- 
sity matrix element pio{t) is given approximately by the 
undamped value 



Using the bootstrap method described in Section Fill Al 
we solve Eqs. (|122|) for the element pio{t) with the initial 
conditions pio(O) = i and p2o(0) = 0. If we again assume ^.^^^ ^ ^^/^^ ^ 2g7r/i.a, 

that the interaction is unsaturated (i.e., |ria|/72o *C 1) 
so that |yOio(i)| ^ |P2o(i)| for all t, and that — to ze- 
roth order in \ila\ — |pio(i)| varies slowly compared to 
+ 720' then Eqs. (jl22|l yields the approximate solu- 



■ 2g7r. 



(127) 



(128) 



tion 



(-710 + -iWio) t 



(124) 



where 



10 



«720 



-7I0 



(125) 



It is straightforward to extend Eq. (|124|l to include a 
coherent unsaturated interaction with N independent 
atoms localized within a volume that is small com- 
pared to (ttw^) /2A at r = 0. We begin by extend- 
ing both the iV-atom density matrix and Hamiltonian 
given by Eq. ifSU]! and Eq. l(5T|) . respectively, to include 
the coherent superposition with the empty resonator, as 
was done in Eq. H113(l and Eq. (|114|l . We find that 
Eq. (|124|l remains unchanged if we make the substitution 
\^a^ N |fia|^ and 710 A^72i/4, so that — in the low- 
dephasing limit, after the transient terms in Eq. H124() 
have decayed — the effect of the placement of N atoms 
in the interaction region is to replace the time t with 
Nt. In this hmit, Eq. (|124|l and Eq. H125|) are entirely 
consistent with Eq. lf55|l and Eq. (|^ . except for the ap- 
pearance of the photon number Ua in the denominator 
of Eq. jnOI)- Since Eq. (|125|l was obtained using an Ua- 
photon Fock state (rather than a coherent state, as was 
implicitly used in the semiclassical estimate of Wq), we 
expect an additional factor of Ua from the analysis lead- 
ing to Eq. l(TTUI) . 

In the unsaturated limit, our analysis of the conditions 
required to obtain a particular phase shift is significantly 
different from that of the undamped case. Clearly, if we 
wish to accumulate a large phase shift before the sys- 
tem state has become significantly mixed, the resonant 
detuning should satisfy the inequality 



« 720 <.ya < J^^\na\ 

V 710 



(126) 



This constraint can only be met in the weak-field case if 
the dephasing rate between the upper and lower atomic 



At time t ^ tq, Eq. (|124|l gives for the magnitude of 
Pio {tq) 



\Pwitq)\ 



■ cxp 



720 \^a\ 



■2q7r 



(129) 



Now, in order to achieve a phase shift of — tt, we must 
choose a long delay time such that q w i^a/2 |ila|^, giving 
tq TTiya/ \^a\ and 



1 



Pwitq) ~ -77 cxp TT 



720 



(130) 



It is clear that we must detune the laser field such 
that Va ^ 720 so that we can minimize the effects of 
decoherence, but it is not obvious how to choose a specific 
value of Va- First, we can define the fidelity (a measure of 
distance between quantum states) of two density matrices 
pi and p2 as0 



F{pi,p2) =Tryp/ P2P1 



1/2 1/2 



(131) 



or, in the case of a pure state and an arbitrary state 
P, 



i^(V,p) = V(^IpIV'). 



(132) 



Applying Eq. H132I) to the density matrix given by 
Eq. H113|) and the pure state 

\4,{tq)) = -^(|l,0,0,na) - |l,na,0,0)). 



we obtain 
1 



1 -|- exp 



720 



720 TT 

'^2' 



(133) 



Second, we can compute the entropy of p{tq) using the 
definition'^ 



^ ^ - Tr [plog2(p)] = - E log2 (A,) , 



(134) 
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where the sum is carried over the nonzero eigenvalues of 
p. Now, the density matrix elements p2i itq) and p20 {tq) 
are proportional to |ila| /i/q, while P22 {tq) is proportional 
to \na\ /vl. Therefore, by Eq. we ignore these 

terms in the density matrix given by Eq. H113|l . and we 
apply Eq. H134(l in the limit 720 / Va ^ 1 to obtain 



S{tq) 



720 



1 



720 



log ( TT 



log2(e) 



(135) 



In principle, we can choose the value of Va to obtain par- 
ticular values of the entropy and fidelity, and then — in 
an TV-atom system — allow the system to evolve until 
time t — tq/N to accumulate a — tt phase shift. In prac- 
tice, in many cases the probe field frequency cannot be 
modified post hoc, particularly in applications where a 
phase shift other than — tt is required and/or more than 
one type of atom or molecule is placed in consecutive 
interaction regions. 

It is already clear from Eq. H135|l and Eq. (|133|l that 
we must have 1^0/720 > 30 if we wish to hold F > 0.95 
and S < 0.5 (base 2) for quantum information purposes. 
If, as an example, we also have \^la \ /720 = 0.2, then after 
a time Nj2Qtq = 7757r we will obtain a linear phase shift 
of — TT. However, if we require F > 0.9995 and S < 0.01, 
then we must have Vall20 > 4000 and wait until a time 
Nj2Qtq = 100, OOOtt to obtain a linear phase shift of —tt. 
In other words, even though we are using a large detun- 
ing, requirements of small entropy and high fidelity imply 



that we will need either extraordinarily long interaction 
times or many identical noninteracting atoms to achieve 
nontrivial phase shifts. 



B. The Four-Level Atom 

A calculation of the real and imaginary parts of pio (t) 
in the case of the four-level atom proceeds in essentially 
the same fashion as the corresponding calculation for the 
two-level case described in the previous section. Again we 
further extend the product state basis given by Eq. 
to include a "second rail" as an alternative quantum path 
for the Ua probe photons, corresponding to a 6 x 6 den- 
sity matrix and Hamiltonian (extended from Eq. H87() 
as Eq. (|114|l was from Eq. (|29|l ). The matrix elements 
Pio(i), P2o(i), /03o(i), and pio(t) are mutually coupled, 
and we seek an approximate solution using the bootstrap 
method used in previous sections. We neglect transient 
(homogeneous) solutions to the coupled ODEs, and we 
assume that pio{t) = ^ is much larger than the magni- 
tudes of the other three elements. Under these condi- 
tions, in the unsaturated limit \il.a\^ /720 <C 1 we obtain 
the quasi-steady-state solution 



(136) 



where 
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710 and 720 are given by Eq. H123|) . and 



where 



730 = 73' and 
1 , 

740 = 2'^4+74- 



(138a) 
(138b) 



If we set = = to minimize absorption, and if the 
dephasing rates are much smaller than the depopulation 
rates of the upper atomic levels, then Eq. (|137|l becomes 
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Note that Eq. (|139bl) has precisely the same form as 
the corresponding two- level result given by Eq. (|125|l . 
Hence, the nonlinear phase shift derived from Eq. E 



can be as large as the corresponding linear phase shift 
obtained from Eq. H125() . indicating the presence of an 
enormous third-order nonlinearity that couples the three 
fields. 3, 53, 62, 64J In principle, this nonlinearity requires 
only modest detunings to provide both a high differential 
phase shift and a low absorption rate in the semiclassical 
realm. 

However, in the quantum information processing ap- 
plications described above, we must also check that the 
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photon-atom system is effectively disentangled when the 
dipole interaction is switched off. In the limit ^ 7201 
the earliest elapsed time ti required to obtain a phase 
shift of — TT and the corresponding fidelity and entropy 
are respectively given by 



F 
S 



par 

and 



720 



log ( ^TT 



log2(e). 



(141) 
(142) 
(143) 



Given a sufficiently long interaction time, it is clear 
from Eq. H14Ua|l that a large phase shift can be accumu- 
lated using a relatively small net detuning i>c even when 
flc represents the Hamiltonian matrix element describ- 
ing the dipole interaction of a single photon with a single 
atomJOT*! However, the value of Vc needed to maintain a 
high fidelity and a low entropy depends on other system 
parameters. For example, in the limit |53|| 
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(144) 



we note that we must have Vc/l20 ~ i^c/lio S> 1 to 
achieve F — > 1 and 5* — > 0. Hence, relatively large 
detunings are still required when the four-level system 
is used for quantum information processing applications. 
However, in the limit where the spontaneous emission 
rate of atomic state |4) in Fig. has been strongly sup- 
pressed. 



in. 
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(145) 



we find from Eq. ifH^ and Eq. that i>c/720 ~ 

i>c/720; a constraint that can be satisfied easily for small 
detunings whenever |rif,|^ / jild^ ^ 1. Suppression of 
spontaneous emission from level |4) can be achieved in 
at least two different ways. First, the interaction region 
can be placed within a photonic bandgap crystal (PBC) 
designed so that photons with frequency uJc must be in- 
jected through a defect in the crystal structure. Second, 
a different atomic system could be chosen with an energy 
level structure similar to that shown in Fig.^J where the 
final dipole transition in Fig. 



7 a 



has been replaced by 
a two-photon transition to the metastable atomic level 
|4). Although the details of the calculations leading to 
Eq. Ijl39b|) will certainly change for this system, the rel- 
ative weakness of the two-photon transition amplitude 
can be effectively offset by a suitably smaller choice of 
the value of the detuning frequency Vc- 

If the controlled couphng transition |2) < — > |3) is 
driven by a Fock state, then the four-level atom + field 
system provides a large cross-Kerr nonlinearity of the 
form -ffKcrr = "-a fic whcn the constraint given by 
Eq. (|145|l is satisfied. Neglecting the imaginary part of 




FIG. 12: Modified interaction between a four-level A/" atom 
and a three-freque ncy e lectromagnetic field. The final dipole 
transition in Fig. |7(a)| has been replaced by a two-photon 
transition to the metastable atomic level |4). 



Eqs. H139|l . we find that the evolution corresponds to that 
of the Kerr Hamiltonian given by Eq. (fTUT|l . with the 
nonlinear coefficient given by Eq. (|105|) . In this case, the 
entanglement of the atoms and fields is negligible, and 
the fidelity of the final state is high. However, if the 
coupling transition is driven by a coherent state, then — 
as shown in Fig. ^| — the electromagnetic intensity of 
that state (i.e., the magnitude of the Poynting vector) 
must be quite large to ensure that the atoms and fields 
are completely disentangled at the conclusion of the gate 
operation. This condition requires that \ai,\^ ^ 1, reduc- 
ing the magnitude of the effective cross-Kerr nonlinearity 
given by Eq. l(TU5|) . 

We have demonstrated so far that it is possible to apply 
an arbitrary phase shift to an initial state cq |0) -I- ci |1) 
of a photonic qubit tuned to the |1) < — > |2) transition 
in Fig. [7| resulting in the final state cq |0) -I- cie"^ |1). 
Clearly it is straightforward to perform Hadamard gates 
on qubits encoded in photons in this manner, through use 
of simple linear optics (beamsplitters). With the tune- 
able EIT phase shift gate, the required range of single 
qubit gates needed for universal quantum information 
processing is therefore covered. The other necessary in- 
gredient for universal quantum processing 0, |3 163, IzOl 
is a two qubit entangling gate. For photonic qubits which 
generally interact very weakly with each other, this is 
the more difficult gate to realize. One solution is to use 
measurement and feedback to create an effective strong 
non-linearity between photonic qubits. |^ However, it is 
clearly of real significance for photonic quantum informa- 
tion processing to consider the possibility of a direct non- 
linear coupling between photonic qubits using an EIT 
system, to realize, for example, a conditional two-qubit 
phase gate. 

Consider a case (depicted schematically in Fig. [T^ 
with two photon number encoded qubits (target and con- 
trol), where the target qubit is tuned to the |1) < — > |2) 
transition frequency LOa and the control qubit is tuned to 



23 



Control 



Control 




Target 



Target 



FIG. 13: Schematic diagram of a two qubit conditional phase 
shift which transforms the state coo |00) + coi |01) + Cio 1 10) + 
cii 111) to coo |00) + coi jOl) + cio |10) + ciie*'^* |11). When 
Xt = TT the two qubit phase gate is implemented. 



the 1 3) < — > 1 4) transition frequency Wc of Fig. [7| so this 
mode is now a quantum rather than classical control field. 
As shown above, if no photon is present in the |3) < — > |4) 
transition, then the target qubit Co |0) + Ci |1) acquires 
no phase shift. However, if a photon with frequency lOc 
is present in the |3) < — > |4) transition, then the target 
qubit evolves to cq |0) + cie*"^ |1). Hence, this system im- 
plements a conditional phase shift and is extremely useful 
for quantum information processing. In fact, for a condi- 
tional phase shift of — tt the operation provides a univer- 
sal two-qubit gate capable of maximally entangling two 
initially unentangled photonic qubits. The input product 
state i(|00)-l-|01)-l-|10)-l-|ll)) can be transformed to the 
maximally entangled state |(|00) + |01) + |10) - |11)), as 
shown schematically in Fig. 1131 Therefore, in principle, 
universal optical quantum information processing can be 
performed with such FIT systems. 

To illustrate this we consider the error probability in 
generating the conditional phase shift on the |11) ampli- 
tude. Fig. shows the error probability in generating 
the phase-shifted amplitude for single photons in the a 
and c modes of the four-level system, as a function of 
the average photon number in the coherent drive applied 
to mode b, for various values of the chosen phase shift. 
Clearly as long as the applied drive in mode b is pushed 
towards being a classical field and so the KIT system pro- 
vides an accurate cross-Kerr nonlinearity, the |11) ampli- 
tude can in principle receive a large and accurate phase 
shift. The dependence of the nonlinearity on the various 
parameters is given in Eq. H108I) . From Fig. we see 
that the gate requires a drive field with a large jabp, but 
from Fq. (|108|l this reduces the strength of the nonlin- 
earity, which increases the value of Nt required for the 
chosen phase shift (for a given detuning). So, as discussed 
in Section IIV Al and Section IIVBI it is crucial that the 
effects of decoherence are kept small, in particular the 
spontaneous emission from level |4), in order to perform 
a two-qubit gate with small error. Detailed analysis of 
the effects of decoherence on the FIT phase gates will be 
addressed in future work. 

Clearly the conditional phase shift can be put to good 
use in other quantum processing applications. One ex- 
tremely useful device is a high efficiency non-demolition 
detector for photons. If instead of a single photon state 



I 



10 



10' 



10 



10-^ 



10 



10 



—1 1 — I I I I 1 1 1 1 1 — I I I I 1 1 



n„=l 
n =1 



-■ i]) = 7t/4 
-■^ = K/2 

■ (]) = 3jt/4 
— (]) = 7t 



_i I ' I I I 



10 



10" 



10 



10 



10 



FIG. 14; Numerical evaluation of the error probability 1 — 
for several values of the net phase shift (j>, assuming that the 
input Fock state has Ua — 1 and ric = 1. We have used 
Eq. llT^ to compare an ideal phase-shifted state Eq. HUGH 
to the actual final state Eq. Hl()7|l . 



on the |3) < — |4) transition, a weak coherent state is 
input, a measurable (by standard techniques) phase shift 
arises conditional on the single-photon Fock amplitude of 
a qubit in the |1) < — > |2) transition. In effect a projec- 
tive measurement in the computational basis can be per- 
formed on the photonic qubit, with the qubit being avail- 
able for re-use afterwards rather than being absorbed. 
This can be achieved with high efficiency (> 99%) and 
with just a few hundred atoms in the KIT system. Details 
of this detector are reported elsewhere. \T\\ 



V. CONCLUSIONS 

In Section IIII CI we demonstrated that it should be 
possible to use the four-level atomic system of Fig. [7| 
as the foundation for both a (so-called) quantum switch 
and a phase-shifter. For example, the quantum switch 
is, in principle, simple to implement: given the availabil- 
ity of a photon with frequency ujb — ^32, if a photon 
with frequency lOc = ^'-'43 is present (i.e., fie 7^ 0), then 
a probe photon with frequency Wa = W21 will be ab- 
sorbed; otherwise, if Q,c = 0, then the probe photon will 
be transmitted. All this shows that conventional classical 
information processing operations are possible on optical 
data — single bit phase shifts and conditional (two-bit) 
switching. These effects have significant potential for ap- 
plications to conventional optical communications. 

In Section Hvl we demonstrated that it should be pos- 
sible to operate the four-level atomic system as a "dual 
rail" photon qubit phase shifter, provided that the spon- 
taneous emission from the atomic level |4) can be sup- 
pressed. The size of the phase shift and the fidelity of 
the gate are quantified in terms of the system parame- 
ters in our model, so the trade-off between the size of 
the phase shift, the accuracy of the gate, the time of op- 
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eration and the atomic and control parameters can be 
analyzed in detail. A single photon phase gate, coupled 
with others that can be performed using linear optical el- 
ements, enables the performance of arbitrary single qubit 
operations. Universal quantum information processing 
requires the addition of a suitable entangling two-qubit 
gate.0, We have also demonstrated that in 

principle, the phase shifter arrangement can be turned 
into such a gate — a two-qubit conditional phase gate — 
by using a qubit input also on the control field at fre- 
quency Wc- Therefore such EIT systems can in principle 
be used to enable universal quantum information process- 
ing with "dual rail" photon qubits. Coupled with ideas 
such as quantum memory for photons |2^ 0; HE US 
and non-absorbing photon detectors, [tJ it is clear that 
EIT systems present a very promising route forward for 
few-qubit quantum information processing. 

Indeed, given the experimental progress with EIT phe- 



nomena over the last few years, it seems likely that these 
QIP applications can tested and developed over the next 
few years. Of course, more detailed research still needs to 
be done. For example, further refinements to our model 
to include coherent wavepacket representations of the 
fields are needed to realistically assess in detail the per- 
formance of the two-qubit gate. This will be addressed 
in future work. 
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